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ABSTRACT 

We construct three extreme different scenarios of the star formation histories apphcable to a 
sample of dwarf galaxies, based either on their present metallicity or their luminosity. The three 
possible scenarios imply different mechanical energy input rates and these we compare with the 
theoretical lower limits established for the ejection of processed matter out of dwarf galaxies. The 
comparison strongly points at the existence of extended gaseous haloes in these galaxies, acting 
as the barrier that allows galaxies to retain their metals and enhance their abundance. At the 
same time our findings strongly point at a continuous star-forming process, rather than to coeval 
bursts, as the main contributors to the overall metallicity in our galaxy sample. 

Subject headings: ISM: bubbles, ISM: abundances, general - starburst galaxies 



1. Introduction 

Given the extreme violence of stellar ejection 
processes, either through supernova explosions or 
strong stellar winds, the presence of strong re- 
verse shocks assures that upon thermalization the 
ejected matter would reach temperatures {T ^ 
10^ — 10^ K) that would strongly inhibit recombi- 
nation, and thus the detection of the newly pro- 
cessed material at optical and radio frequencies. 
Furthermore, it is now well understood that this 
hot high pressure gas fills the interior of superbub- 
bles and drives the outer shock that sweeps and 
accelerates the surrounding ISM into an expand- 
ing large-scale supershell. The continuous energy 
input rate from coeval starbursts or continuous 
star forming phases, due to last until the last 8 
M0 star explodes as supernova, reassures that the 
high temperature of the ejected matter is main- 
tained above the recombination limit (T ^ 10^ K) 
allowing superbubbles to reach dimensions in ex- 
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cess of 1 kpc. During this time the highly metal 
enriched superbubble gas does not diffuse either in 
the expanding shell nor in its immediate surround- 
ings. It has been convincingly argued in recent 
years that the metallicity detected in star-forming 
galaxies results from their previous history of star 
formation and has nothing to do with the metals 
presently ejected by their powerful starbursts (see 
Tenorio-Tagle 1996). One in fact would have to 
wait until the end of the supernova phase (> 40 
Myr) for the full ejection of all heavy elements syn- 
thesized by the massive stars. At this stage up to 
almost 40% of the starburst mass would have been 
returned to the ISM and heavy elements, such as 
oxygen, would have reached their expected yield 
value. 

Note that the newly processed heavy elements 
are generated within a region of about 100 pc (the 
typical size of starburst) while upon cooling of the 
superbubble interior and their fall back towards 
the galaxy, the heavy elements are disseminated 
over a region of several kpc in size. After such 
a large-scale dispersal, diffusion leads to a thor- 
ough mixing of the heavy elements within the host 
galaxy ISM, enhancing its metallicity. 

Indeed the final step towards mixing within 
the ISM is promoted by diffusion since diffusion 
is largely enhanced in high temperature gases. 



Therefore, if a new major centre of star forma- 
tion develops and causes the formation of an HII 
region, the heavy elements will rapidly mix to pro- 
duce a uniform abundance. In the absence of star 
formation however, the heavy element droplets 
will also diffuse but over larger time-scales (see 
Tenorio-Taglc 1996). All together the process of 
mixing takes several 4-6 x 10^ yr to promote 
an enhanced and almost uniform abundance in a 
region of up to several kpc in radius around the 
original star burst. 

On the other hand, the observational evidence 
of powerful starbursts in dwarf galaxies has led to 
the idea that, due to their rather shallow gravita- 
tional potential, supernova products and even the 
whole of the interstellar medium, may be easily 
ejected from the host dwarf systems, causing the 
contamination of the intra-cluster medium (Dekel 
& Silk 1986; De Young & Heckman 1994). How- 
ever, the indisputable presence of metals (in what- 
ever abundance) in galaxies implies that the su- 
pernova products cannot be completely lost in all 
cases. 

This issue has been recently addressed by Silich 
& Tenorio-Tagle (2001; hereafter referred to as Pa- 
per I) and D'Ercole & Brighenti (1999) who show 
that within a dark matter (DM) and a low-density 
extended gaseous halo, a galactic wind does not 
easily develop, making the removal of the galaxy's 
ISM very inefficient. Here we compare the theo- 
retical estimates with some well known starbursts 
in blue compact dwarf (BCD) galaxies. In partic- 
ular we derive three different possible star forma- 
tion history scenarios that assume either a very 
young coeval starburst or an extended phase of 
star formation of 40 Myr or 14 Gyr, to account for 
the inferred star formation rate (in the first two 
cases) and the observed mctallicity (in the third 
one) from our sample galaxies. For these galaxies 
we know: their HI mass, metallicity and Ha or 
H/3 luminosities, and these allow us to infer the 
energy input rate expected in each of our galaxies, 
for each of the three likely histories of star for- 
mation (see Section 2). Section 3 is devoted to a 
comparison of the energy input rates derived for 
our galaxies with the recent estimates of the rates 
required to expel either metals or the ISM from 
galaxies containing a total ISM mass in the range 
of 10*^ to 10^° Mq (see Paper I). Section 4 gives a 
discussion of our results. 



2. The Star forming history of BCD galcix- 

ies 

Our sample of dwarf galaxies (sec Table I) was 
selected from Papaderos et al (1996), Mas- Hesse 
& Kunth (1999), Marlowe, Meurer & Heckman 
(1999), Lequeux et al. (1979). We took care to 
exclude galaxies which could be tidal dwarfs, be- 
cause in such a case their metallicity would not 
be the result of the evolution of the present day 
system. The sample also includes a few dirr galax- 
ies from Skillman, Kennicutt & Hodge (1989). In 
Table 1 the first column provides the galaxy iden- 
tification and columns 2 and 4 give the HI mass 
and metallicity as reported in the references given 
in columns 3 and 5, respectively. The fact that 
many galaxies, such as IZw 18 (Vilchez & Iglesias- 
Paramo, 1998; Legrand et al. 2000; see also van 
Zee et al. 1998), IIZw 40, NGC 4214 (Kobulnicky 
& Skillman 1996), NGC 1569 (Devest et al. , 
1997), and others, present a uniform metal abun- 
dance within their HII region gas, irrespectively of 
the chosen slit location, and the arguments given 
in section I regarding the unlikely possibility that 
the metals come from the present generation of 
stars, have led us to assume that galaxies have 
been evenly polluted throughout their histories, 
up the levels now displayed in their HII region gas. 

We examine three different scenarios of the star 
formation using the observed properties of our 
sample of dwarf galaxies. For our first two mod- 
els, the observed H^, luminosities (column 6 of Ta- 
ble 1) were used to estimate the gas mass con- 
verted in stars. We derive the UV photon output 
from the exciting cluster using Leitherer & Heck- 
mann (1995) (hereafter referred to as LH95) rela- 
tion L{Ha) erg s^^ = 1.36x10-'^^ N{H'^) s^^, and 
then convert it into the SFR. References to the 
observed Ha luminosity are given in column 7 of 
Table 1. These correspond to the total luminosity 
obtained from digital photometry on CDD images 
using an Ha filter. Whenever the Ha luminosity of 
a galaxy was not available, we used H/3 instead and 
using the relation: Ha/H/3 = 2.857, as in LH95. 
Galaxies for which the Ha luminosity has been in- 
ferred from their H/3 luminosity are marked with 
an asterisk in the column 7 of Table 1. In these 
cases, it can be assumed that the observed H/3 flux 
also corresponds to the total flux, either because 
the slit apperture has been chosen large enough 
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for a realistic comparison with lUE data (Mas- 
Hesse & Kunth, 1999, Storchi-Bergmann, Kinney, 
Challis, 1995), or photometrically derived from 
CDD images (Vilchez, 1995 Gallagher, Hunter & 
Bushouse, 1989). 

We have first assumed that all stars have re- 
sulted from a coeval burst. Under these condi- 
tions, it is well known that the mechanical input 
rate through winds and supernovae leads to an al- 
most immediate constant value to be maintained 
throughout the supernova phase (4 xlO'' yr). On 
the other hand the ionizing flux is to remain con- 
stant at its maximum value only for the first 3.5 
Myr, and then decreases as (see Beltrametti 
et al. 1981). The implication of this rapid de- 
cay is that after lO'' yr the UV photon output 
is two orders of magnitude smaller than its origi- 
nal value, and thus the HH region phase is much 
shorter than the supernova phase. Here we assume 
that our clusters are indeed coeval and young (< 
3.5 Myr). The Ha luminosity is also in this case 
a direct indicator of the number of ionizing pho- 
tons. However, to infer the mass of the starbursts, 
the number of ionizing photons should be linearly 
scaled using Figure 37 of LH95, which gives a value 
of N{H°) ~ 10^3 photons s"! for every 10^ Mq of 
stars formed. The corresponding mechanical en- 
ergy input rate can thus be evaluated using the 
calibration given by these authors for an instan- 
taneous burst of 10^ Mq (LH95, their Figure 55). 
The values found for the mass in stars and me- 
chanical energy injection rate are given in column 
2 and 3 of Table 2, for each entry. 

In this model, our mass estimates will be wrong 
by a factor of ten if the age of our starbursts 
would be 6.5 Myr, and by a factor of 100 if 
they all are 10 Myr old. The presence of strong 
Wolf-Rayet features in the spectra of many of the 
sources in the sample, as in IZw 18, (see Kunth & 
Ostlin 2000, and references therein) certifies how- 
ever their young age (^3.5 Myr). 

In our second model the present star formation 
rate was assumed to be continuous for the last tens 
of Myr. In this case, the ionizing photon output 
reaches a maximum constant value after 3-4 Myr 
while the mechanical luminosity of the starburst 
rises to reach its asymptotic maximum value only 
after 40 Myr (LH95). We thus use their continuous 
star formation model to derive the SFR from the 
number of ionizing photons {N{H^)) as well as the 



corresponding energy injection rate, assuming a 
continuous and constant SFR over the last 40 Myr 
(LH95, Figures 38 and 56). The resultant SFR and 
the corresponding mechanical energy input rate 
for each of our galaxies are given in columns 4 and 
5 of Table 2. Note that despite the fact that there 
is no information on the time since star formation 
began in our galaxies, our assumption that t > 40 
Myr, leads to an upper limit for the mechanical 
energy input rate as implied from the derived SFR. 

In models 1 and 2, we have assumed that all 
the ionizing photons produced by the stellar clus- 
ters are completely used up in the ionization of the 
HII region gas. Note however that if a fraction of 
these, say 50%, manage to freely escape the neb- 
ulae, or are absorbed by dust, our final estimates 
are only off by a factor of two. 

For our third model we have evaluated the 
amount of stars required to produce the observed 
oxygen abundance in each galaxy of our sam- 
ple. Because this approach is based on the entire 
galaxy evolution, we used here a closed box model 
and the instantaneous recycling approximation de- 
spite of the time delay between the metal ejection 
and their complete mixing with the galaxy ISM. 
In the low metallicity limit (or for large remaining 

gas fraction u = ^1°.°^ ) the ISM gas metal- 

licity Z is given by: 

z = -y (1) 

where y is the net yield (Maeder, 1992). This ap- 
proximation is valid for all the objects in our sam- 
ple. For oxygen this relation becomes: 

O _ 
H + He + metals 

where yo is the oxygen yield. As we only have a 
measurement of the hydrogen content, we have as- 
sumed a constant He/H ratio: Mgas = Mhi+Hb = 
1.34 Mhj (Kunth & Sargent 1986). 

Thus, reverting equation 2 we can derive the 
mass of stars required to produce the observed 
oxygen abundance in our objects, as a function 
of their HI mass and the oxygen yields. 

It is important to notice that the derived mass 
in stars is a lower limit. Indeed, under the as- 
sumption of a closed box model, all of the metals 
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produced remain within a galaxy and contribute 
to globaly enhance its metalhcity. If some of the 
metals could escape into the intergalactic space, 
the amount of stars needed to produce the actual 
measured metallicities will have to be larger. Of 
course, this mass also depends on the yields used, 
and the normalization of the yields depends on the 
adopted initial mass function. Here we used the 
oxygen yield Yq = 0.0026 computed from Maeder 
(1992) models for an initial mass function with a 
power law index -2.35 between 0.1 Mq and 120 
Mq. The evaluated stellar mass in the range 1- 
120 Mq needed to produce the observed global 
metallicity of our galaxies is given in column 6 of 
Table 2. 

In the last scenario we have assumed that star 
formation has taken place all over the volume of 
the host galaxies. We have then divided the de- 
rived stellar masses by 14 Gyrs to compute the 
continuous star formation rate required to produce 
the observed metallicity over a Hubble time. This 
star formation rate is given in the column 7 of Ta- 
ble 2. Note that this scenario leads to a lower limit 
for the mechanical energy injection rate. Indeed, 
if the same amount of stars were formed during 
burst events, the expected mechanical energy in- 
put rate should become momentarily larger. 

The corresponding mechanical energy injection 
rate (column 8 of Table 2) was computed using the 
normalization log{L„iech) erg = 41.8. This 
is a mean value predicted by a continuous star- 
burst model, for different star metallicities, after 
the equilibrium regime is reached (see Figure 56 of 
LH95) under the assumption of a Salpeter initial 
mass function between 1 and 100 Mq and star for- 
mation rate of 1 Mq yr^^. Therefore we have used 
this value for our sample regardless of the galaxy 
metallicity. 

Note that if a continuous star formation pro- 
cess, spread over the whole galaxy volume and 
over a Hubble time, is required to produce the 
observed global abundances detected in our galax- 
ies, then the implied SFRs and total mechanical 
luminosities are comparable to the values derived 
for either our model 1 or 2. This is remarquable 
since for models 1 and 2, the SFR, the mass in 
stars, and the mechanical luminosities were in- 
ferred independently from the present-day galaxy 
luminosity. Although metals presently being pro- 
duced are hidden from view (at least in the optical 



regime) because the recombination process is in- 
hibited at high temperatures when crossing the re- 
verse shocks (see section I), we can derive, for each 
galaxy, the number of "typical" starburst events 
similar to the present ones, required to account 
for the observed metallicities. If one divides the 
stellar mass implied by chemical evolution (col- 
umn 6 in Table 2) over the starburst mass derived 
in scenarios 1 and 2, one would find the number of 
events needed to explain the metal content, evenly 
spread over the galactic volumes. Column 9 of Ta- 
ble 2 provides the number of extended (40Myr) or 
coeval bursts, required to account for the metal- 
licity of the galaxies in our sample. Note that the 
derived numbers of bursts are also a lower limit 
estimate as the closed box model provides a lower 
limit on the stellar mass. Nethertheless, the large 
number of bursts required in most of the cases, 
both in models 1 and 2, simply reflects the fact 
that it is the whole galactic volume what aught to 
be contaminated. These large numbers seem un- 
realistic and thus one is lead to conclude that the 
main agent causing the global metallicity present 
in galaxies is a uniform, constant and long-lasting 
SFR, spread (as in model 3) over the galaxy vol- 
ume. Starbursts may take place throughout the 
history of the galaxy, with a star-forming activ- 
ity similar to what is found in the whole galactic 
volume, however, their contribution to the total 
metallicity would not be the most significant. 

Finally, we compared the values of the derived 
mechanical energy injection rate in the three con- 
sidered cases, with the predictions of the hydrody- 
namical models of McLow & Fcrrara (1999) and 
Silich & Tenorio-Tagle (2001). The results are 
shown in Figure 1. 

3. Comparison with the theoretical esti- 
mates 

3.0.1. The energy requirements 

Figure 1 shows our energy estimates resultant 
from the numerical integration of the hydrody- 
namic equations for a coeval starburst with a con- 
stant energy input rate during the first 4 x 10^ yr 
of the evolution, acting on an ISM density distri- 
bution with a total mass in the range of 10^ to 10^^ 
Mq. Two extreme density distributions were as- 
sumed: either rotating flattened disks (lower two 
lines) or spherical galaxies without rotation (up- 
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per two lines; sec Paper I). Clearly disk- like mod- 
els require less energy to eject their metals into the 
intergalactic medium. This is because the amount 
of interstellar gas that has to be blown out of a 
disk-like galaxy to open a channel for metal ejec- 
tion into the intergalactic medium is much smaller 
than in the spherically symmetric limit, where all 
of the metal-enriched ISM has to be accelerated to 
reach the galaxy boundary. 

In Figure 1, values below each line imply total 
retention, while the region above each line indi- 
cates the expulsion of the hot superbubble interior 
gas (with the new metals) out of disk-like distri- 
butions, and of the new metals and the whole of 
the ISM in the spherical cases. Each line that 
separates the two regions marks the minimum en- 
ergy input rate needed from a coeval starburst to 
reach the outer boundary of a given galaxy with a 
supersonic speed, regardless of the time that the 
remnant may require to do so (see Paper I). 

The energy input rates derived for our sample 
of galaxies (see Table 2) in each of the three star 
formation scenarios here considered have been in- 
corporated in panels la - Ic. Note that the energy 
input rate derived in model 3 from the continu- 
ous star formation activity able to account for the 
metallicity of the host galaxies, although originally 
assumed to take place all over the galaxy volume is 
here compared with the theory estimates derived 
for spatially concentrated starbursts. 

4. Discussion 

Based on sound observed facts such as the HI 
mass, metallicity and either the Ha or H/3 lumi- 
nosity of star-forming dwarf galaxies we have con- 
structed three different likely scenarios for their 
star formation histories. The various extreme pos- 
sibilities result from considering either young co- 
eval starbursts or extended phases of star forma- 
tion lasting 40 Myr and 14 Gyr respectively. 

For the coeval starburst scenario, as well as for 
the second extended and constant star formation 
rate model, the star formation rate and the to- 
tal mass of the starburst have been inferred from 
the observed Ha and H/3 luminosity of each of the 
galaxies. The third scenario assumes a constant 
star formation rate over a Hubble time (14 Gyr), 
as low as it is required to account for the observed 
metallicity in each of these galaxies. Based on the 



synthetic properties of coeval and extended star- 
bursts we have inferred the energy input rate to 
be expected in each of the cases and these have 
been compared with the recent theoretical esti- 
mates of the minimum energy input rate required 
to vent the processed matter into the intergalactic 
medium (see Paper I). 

In all three star formation history scenarios, all 
galaxies in our sample lie above the lower limit first 
derived by Mac Low & Ferrara (1999), and then 
confirmed in Paper I, for the ejection of metals 
out of flattened disk-like ISM density distributions 
when energized by relatively low mass starbursts. 
Under these assumptions, the observed metallicity 
in these systems cannot be accounted for. Our re- 
sults however, place the bulk of our sample, in all 
three considered star formation history scenarios, 
below the minimum energy input rate range of val- 
ues derived for galaxies that present an extended 
low density gaseous halo. And thus, as shown in 
Paper I, it seems to be the mass of the extended 
low density halo, the one that acts as the barrier 
against metal ejection into the IGM. 

Note that in the 14 Gyr extended star forma- 
tion case, all galaxies lie below the lower limit im- 
posed by the presence of an extended low density 
halo (see Figure Ic). A few objects namely: NGC 
1569, IZw 18, Haro 2, and Mrk 49 surpase this 
lower limit in the case of a 40 Myr extended star 
formation history (see Figure lb), as it is also the 
case of NGC 1569 in the coeval starburst scenario 
(see Figure la). These facts imply that the limit 
derived by Silich & Tenorio-Tagle (2001) is indeed 
a lower limit as it does not consider all possible 
phases of the ISM. nor magnetic fields and it has 
(as other models in the literature) assumed a con- 
stant energy input rate. All of these are issues 
that if thoroughly considered could substantially 
raise the energy input rate required for mass ejec- 
tion out of dwarf galaxies hence the metals are 
indeed trapped in all cases. A second possibility 
is that if the metallicity in these galaxies results 
from several extended star formation phases, or 
from several coeval bursts (sec column 9 of Table 
2), the metals produced by the present phase of 
star formation will only in a few of the cases be 
able to escape into the IGM. 

However, given the large number of star for- 
mation events that have been derived for many of 
our galaxies in both models 1 and 2, it is obvious 
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that the main contribution to the global metal- 
licity measured in star-forming dwarf galaxies is 
due to a long-lasting and continuous star form- 
ing process, spread over the whole galaxy volume. 
A few coeval and massive bursts, similar to the 
ones presently energizing the galaxies, may oc- 
cur throughout history. Taking 100 Myr time in- 
terval between the sequential starbursts, one can 
find, that in most of the galaxies from oursample, 
instantaneous starburst may account for 1/10 or 
even much less of the detected metal abundances. 

The main implication of our results is that one 
has to invoke the existence of extended gaseous 
haloes in order to account for the metals present in 
galaxies. Predicted haloes, despite acting as a bar- 
rier to the loss of the new metals, have rather low 
densities (< rihaio >~ 10"'^ cm""^) and thus have a 
long recombination time {tree = ^ / {ctnhaio)', where 
a is the recombination coefficient) that can easily 
exceed the lifetime of the HII region {Ihii = 10^ 
yr) developed by the starburst. In such a case, the 
haloes may remain undetected at radio and optical 
frequencies (see Tenorio-Tagle et al. 1999), until 
large volumes are collected into the expanding su- 
pershells. Note that the continuous Vl shape that 
supershells present in a number of galaxies, while 
remaining attached to their central starburst, as 
well as their small expansion velocity (comparable 
or smaller than the escape velocity of their host 
galaxy) imply that the mechanical energy of the 
star cluster is plowing into a continuous as yet un- 
detected medium. Furthermore, in the presence of 
a halo, it is the mass of the halo what sets the lim- 
iting energy input rate required for mass ejection 
into the IGM, and not the mass of the disk-like 
component. This argument applies to all galaxies 
whether spirals, amorphous, irregulars or dwarfs. 
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Table 1: Star- forming dwarf galaxies. Observational data 



Name 


M ( H I ) 
(10'^ 


rcf 


12 + Log(0/H) 


rcf 




L(Ha) 

erg s "'" 


rcf 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




5.4 


1 


8.37 


7 


7 


14 X lO'^" 


9* 


Mkn-49/131S 


1.9 


1 


8.28 


14 


3 


71 X 10^0 


10* 


NGC-2915 


5.9 


IS 


8.4 


7 


1 


30 X lO''^ 


8 


Haro-2 


4.6 


1 


8.4 


6 


1 


20 X lO'^^ 


6* 


II-Zw-70 


3.4 


1 


8.07 


6 


2 


23 X 10^0 


6* 


I-Zw-18 


0.26 


6 


7.18 


6 


4 


29 X 10'^^ 


6* 


Haro-14 


3.2 


1 


8.4 


13 


9 


03 X 10'^^ 


9* 


IC^IO 


13 
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Table 2: Star-forming dwarf galaxies. Starburst properties 
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Fig. 1. — Energy estimates. The log of the critical mechanical luminosity (left-hand side axis) and mass of 
the star cluster MgB (right-hand side axis), required to eject matter from galaxies with a Mjsm in the range 
10^ - 10^° Mq. The lower limit estimates are shown for galaxies with extreme ISM density distributions: 
flattened disks (lower two lines), and spherical galaxies without rotation (upper lines), for two values of the 
intergalactic pressure PiGu/k = 1 cm~^ K (solid lines) and PiGm/k = 100 cm~^ K (dashed lines). Each line 
should bo considered separately (see Paper I) as they divide the parameter space into two distinct regions: 
a region of no mass loss that is found below the line and a region in which blowout and mass ejection occur 
is found above each line, a) - c) Estimates for models 1-3, respectively. 
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